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Abstract 

We consider configurations of D4 — D8 — D8 branes which correspond to large iV QCD 
with non-vanishing temperature and chemical potential for baryon number and isospin. 
We study the holographic dual of this model and find a rich phase structure. The phases, 
distinguished by the values of quark condensates, are separated by the surfaces of first 
order phase transitions. The picture is in many respects similar to the expected phase 
structure of QCD in the chiral limit. 
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1. Introduction and summary 

String theoretic constructions aimed at describing the theory of strong interactions, 
QCD, are becoming more and more elaborate. There are now many versions of the 
gauge/string correspondence where many constraints of the M = 4 Super Yang Mills 
are relaxed, although one still has to consider the limit of large number of colors N. How- 
ever the description of large N QCD-like (asymptotically free) theories proves to be a hard 
task. Typically, the regime where the gauge theory is modified at energies much higher 
than the dynamically generated scale, is not easily tractable on the dual, string theoretic 
side. 

Despite these limitations, holographic models are still useful because very often they 
capture physics qualitatively correct. Hence, they provide a unique window into the dy- 
namics of strongly coupled gauge theories. With more luck, it might be possible to extract 
quantities that are universal, i.e. do not significantly change as one interpolates between 
the holographic dual and QCD. A good place to search for such universal quantities is the 
phase diagram of the theory, especially if the second order phase transitions are present. 

A model proposed in [1] building on [2] which involves a certain configuration of 
DA — D8 — D8 branes is particularly interesting because it gives rise to large N QCD 
with quarks in a certain region of the parameter space. This model, also sometimes called 
"holographic QCD", has a holographic dual which reproduces some expected physics of 
the gauge theory. Chiral symmetry breaking and derivation of the pion lagrangian have 
been described in the original paper by Sakai and Sugimoto [1]. The model has been 
also investigated at finite temperature [3,4]. Nonzero baryon chemical potential has been 
considered in [5,6] where the phase diagram was shown to contain a line of first order phase 
transitions in the \ib — T plane. 1 

In this paper we focus our attention on the case of two flavor (Nf = 2) large N 
holographic QCD with non-vanishing isospin chemical potential, whose value we denote 
by QCD with \ii ^ has been studied in the variety of models (see [13-16] for 
an incomplete list of references.) The phase diagram enjoys rich structure, with first and 
second order phase transition lines separating various phases. This therefore is a promising 

1 In the first versions of [5] the solution with broken chiral symmetry contained an unphysical 
charge. This has been corrected in [6] and in the most recent version of [5]. Recent work which 
discusses related physics and possible additional phases includes [7-11]. The isospin chemical 
potential in the D3 — D7 system was recently discussed in [12]. 
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setting for the search of universal quantities. Besides, theory with non-vanishing pi is 
amenable to lattice simulations and does not suffer from the determinant sign problem, 
which complicates obtaining lattice results for the ^ case. Finally, non- vanishing 
isospin chemical potential can be experimentally relevant. 

The holographic dual of the DA — D8 — D8 branes involves the warped product of 
four-dimensional Minkowski space, a circle parameterized by X 4 e [0, 27ri?4), a holographic 
coordinate U > Uk and a four-sphere (more details are provided below). At non- zero 
temperature the phase where gluons are not confined is described by a black hole metric 
which is thermodynamically preferred for temperatures higher than I/271R4. Most of 
the analysis below will be restricted to this phase. We comment on the appearance of 
confinement/deconfinement transition in the phase diagram in Section 4. 

Fundamental matter is described by Nf pairs of D8 — D8 branes with asymptotic 
separation L. The low energy degrees of freedom which live on the intersections of the D8 
and D8 branes with the DA branes are four-dimensional massless left and right-handed 
quarks. This brane setup thus describes large N QCD in the chiral limit (in the certain 
region of the parameter space). At generic values of baryon and isospin chemical potentials, 
the D8 branes exist in two phases. In the curved phase D8 and D8 branes connect and 
hence chiral symmetry is broken. In the straight phase, the flavor branes are located at 
constant values of X 4 and fall into the horizon of the black hole. This phase of restored 
chiral symmetry is preferred at larger values of T and pb- 

The phase diagram of the system depends on dimensionless parameters LT and pi = 
/J,i/E where E = AttX/L 2 sets the scale of chemical potential. (A is the t'Hooft coupling 
of five-dimensional gauge theory on the DA branes, which becomes four-dimensional after 
Kaluza-Klein reduction.) Fig. 1 contains the picture of the phase diagram in the limit 
LT^O. (It is instructive to compare it with Fig. 1 in [14].) The fundamental domain can 
be taken to be 

-A»i < A*2 < A*i; A*i > (1.1) 
The phases are distinguished by the values of the condensates: 

a, = (V^>, i=l,2; p = (^^rfjipj) (1.2) 

where r 3 is the Pauli matrix in the flavor space. In the upper right triangle defined by 
A*i 5 A*2 > A*c = 0.11 all condensates vanish. Curved blue line separates the phase with non- 
zero pion condensate p at small ps (below the line) from the phases with <j\ 7^ 0, 02 = 
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for —jj, c < p 2 < p c and cri j2 = for p 2 < p c - All these features are also present in Fig. 
1 of ref. [14]. What is not present in our analysis is the second order phase transition of 
[14] which separates the phase with p ^ at small Jij from the phase with cr 1;2 , p = at 
larger pj. This line might be an artifact of the random matrix model used in [14]. 2 

So far we discussed the LT^O case. The diagram does not qualitatively change as 
we increase LT (see Figs. 5,6), although the value of p c decreases until it reaches zero 
at LT ps 0.15 (Fig. 4). This corresponds to the temperature of chiral phase transition at 
Pb = 0. The full phase diagram should be drawn in three dimensions, with Figs. 1,2,5,6 
representing slices at two different values of LT. The surface of the first order phase 
transitions which is represented by the blue curves on the pictures of these slices looks 
in fact somewhat similar to a hyperboloid, with pj axis being the axis of rotation. (Of 
course, there is no rotational symmetry in the pb — T plane in our case) 

The rest of the paper is organized as follows. In the next section we set up the notation, 
and discuss general expectations. We then consider LT-^0 case, establish the necessary 
holographic dictionary and present some analytic and numerical results. In Section 3 we 
consider finite temperature case, where most of the analysis is numerical. Here we sketch 
the form of the three-dimensional phase diagram. We conclude in Section 4. 

2. Nonzero chemical potentials; small temperature 

We consider SU(N) gauge theory with Nf = 2 massless quarks (each described by a 
single Dirac spinor ipi, i = 1, 2) in the large iV limit. Various phases of the theory can be 
distinguished by the values of the condensates (1.2). The matter part of the Lagrangian is 

A = E + iA„) + IH-A^i (2-1) 

i=l,2 

where sum over color indexes of SU(N) is implied. Baryon and isospin chemical potential 
are defined to be 

1 1 

^b = ^i+^2); ^7=2(^1-^2) (2.2) 

respectively. At pi : 2 = the Lagrangian has SU(2) x SU(2) global symmetry. One can 
use the axial part of this symmetry to rotate the cr 12 condensates into p; these condensates 

2 In the view of our motivation, the absence of the second order phase transition is a slight 
disappointment . 
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develop nonzero expectation values, and spontaneously break chiral symmetry. The cor- 
responding Nambu-Goldstone model is the pion. At non-zero value of Hi (and sufficiently 
small ^1,2) the system is supposed to prefer nonzero value for p together with <j\^ = 0. 
(see e.g. [14].) 

This large N QCD has a holographic dual, which is obtained by going to the near- 
horizon geometry of the DA branes on a circle of radius R4 [2] and adding Nf = 2 D8 — D8 
pairs [1]. The model deviates from QCD at energy scales ~ I/R4 which is related to the 
dynamically generated scale as 

1 1 

Aqcd ~ -5-e A 4 (2.3) 

-Tt4 

where A 4 = X/R 4 is the four-dimensional t'Hooft coupling at the cutoff scale. As custom- 
ary in the holographic models, when the scales are widely separated, A4 is small, which 
corresponds to the strong curvature regime in string theory. This is beyond the reach of 
the present string-theoretic technology. We consider instead the opposite regime of large 
A4 where string theory reduces to supergravity + low energy D8-brane dynamics. 

In this Section we restrict to the regime of small temperatures, We assume 

that confinement energy scale is much smaller than the scale associated with the dynamics 
of fundamental matter (I/R4 <C l/£), and even at these small temperature glue is decon- 
fined. 3 (It is not hard to reinstate the effects of confinement/deconfinement transition on 
the phase diagram - we do it in Section 4. ) 

The euclidean metric of the holographic dual is given by 

3 3 

ds 2 = (|) 2 ( Vf ^dXW + (dX 4 ) 2 ) + 2 [(dU) 2 + U 2 dQ 2 } (2.4) 

There is also non-trivial dilaton 

e*=9,(lf (2.5) 
The length scale R is related to the parameters of the gauge theory as 

R 3 = ttA = 7ig s N (2.6) 

Here and in the rest of the paper the string length l s = 1. The dynamics of the fundamental 
matter is described by the DBI action of Nf = 2 D8 — D8 pairs. Asymptotically (as U— >oo) 
the branes and antibranes are separated by coordinate distance L in the X 4 direction. 



Formally, this limit is equivalent to the NJL limit studied in [17-19] 
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They may connect at smaller values of U. This configuration spontaneously breaks chiral 
symmetry. Gauge fields on the D8 branes are holographically dual to the currents of global 
SU(2) x SU(2) symmetry. As in [6] it is convenient to go to the gauge A v = and Wick 
rotate Aq— HAq. Then, the non-zero chemical potential in the lagrangian (2.1) corresponds 
to non-trivial boundary conditions on the values of Aq: 



where Aoi is the abelian part of the field on the i-th brane and antibrane. 

In the absence of external charges, baryon chemical potential cannot affect the curved 
solution [6]. Hence, for the purpose of computing the action, we can set \i\ = —^2 = Hi 
in this phase. There are now two possible curved brane solutions. One solution involves 
the brane with Aq(U '^oo) = ±// connecting with the antibrane with Aq(U ^oo) = ±/i. 
This solution has zero electric field on the resulting curved flavor branes and non-zero 
condensates o\ = 02 • As explained in [17] the holographic dual of the chiral condensate 
involves a non-trivial profile of the open string tachyon and hence is difficult to analyze in 
the DBI approximation. 

Another solution involves brane with Aq(U ^00) = ±/i connecting with the antibrane 
with A (U— >oo) = =F/U. There is nonzero electric field of opposite direction, but same 
magnitude on two curved branes so the branes overlap again. This solution is obtained 
from the previous one by chiral rotation, and hence corresponds to nonzero value of p. 

In the following we will use the fact that in the connected phase (djjAoi) 2 = (djjA 2) 2 
and Xf — X2 where the subscript corresponds to the isospin. The DBI action, up to an 
overall constant, is given by 



where a factor of two in front of the action counts the number of flavors and is added for 
convenience. (Since the only non-vanishing component of the gauge field is A (U), the 
Chern-Simons term vanishes.) Lagrange-Euler equation of motion for the gauge field on 
the branes implies the existence of the conserved quantity, 



A 0i (U^oo) = Hi\ 



% = 1,2 



(2.7) 




(2.8) 



uHduAo) 



= c 



(2.9) 



V / 1-(%A ) 2 + (5) 3 (%XT 
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Eq. (2.9) can also be written as 

,3 

Likewise, the equation of motion for X 4 can be written as 



(^o)W 1 + (W ! ( 2 ,0, 



(8 X^- UjjUl + c 2 )^ r21l v 

[OuJL ) - u6{u5 + _ ^3^5 + c2)t/3 [Z.Li) 

where Uq labels the turning point: it is the value of the U coordinate where djjX^^oo. 
To obtain (2.11) we expressed (<9[/X 4 ) in terms of the constants of the equations of motion 
and then demanded duX^^oc at U = Uq. Substituting (2.11) back into (2.10) gives 

c 2 U 3 

{duAo)2 = U3{US + US<?)-U*(l + <?) (2 - 12) 

where we introduced c via c 2 = UqC 2 . Let us consider the phase with d\jA\ = —duA 2 . We 
need to re-express the values of the integrals of motion through /ii ]2 and L. Starting with 
the former, we have 

\W~V2\= I dU\duA Q \ = U Q F<d (2.13) 

1 JU 

where we used (2.12) in the second equality. We also pass to the rescaled quantities and 
define F(c) as 



x^c 1 



x 3 (x 5 +C 2 ) - (1 + C 2 ) 
Note that F(c) is a monotonic function which satisfies F(0) = and 



(2.14) 



F(Z)~c + o(Z), c<l (2.15) 
In the opposite regime of large c 

F(c) « cict, ci = r , ^ = r ^^^ (2.16) 
Similarly, the separation of the branes can be computed as 

r roo I \ 

-=/ dUduX 4 = \^G{c) (2.17) 

^ JUo V U 

where 



oo 



G ® = k d ^iv^F^ (2 - 18) 



which is a monotonically increasing function with 



G (0) = t^M, J im G(e) = h£l® 
r( 16 ) c^oo r ( 6 j 



(2.19) 



We can use (2.17) to re-express 



U = E G 2 (Z), E o = ^ 



(2.20) 



and then substitute the result into eq. (2.13) 

1 



p 1 -p 2 \ = F{c)G\Z) 



(2.21) 



where p = p/Eq. Eq. (2.21) determines the value of c in terms of p. Next we need to 
compute the action. To regularize, we subtract the value of the action from that of a pair 
of straight branes (duX A = 0) with vanishing chemical potential. 4 The resulting quantity 
is 



5S = 2 



U 



dUU^ + 



J™ (ui-Ui^jl- (duAo)* + (dux*)* 



(2.22) 



With the help of (2.20) this can also be written as 



5S = 2E]>G 7 {c)H(c) 



where 



H(c) = - + J x* 



1 - 



(2.23) 



(2.24) 



x s (x 5 + c 2 ) - (1 + C 2 ] 

This should be compared with the action for the phase where the branes are curved but 
have zero electric field. (This is the phase with cr 1)2 ^ 0, p = 0.) 



5S = 2E*G 7 (0)H(0) 



(2.25) 



The thermodynamically preferred phase has a larger value of SS. One can check that 
(2.23) is larger or equal to (2.25), and hence as soon as pi ^ the phase with p ^ is 



4 Another popular regularization scheme involves introducing an ultraviolet (large U) cut- 
off for the integral. Since we are interested in the difference between the actions for different 
configurations, these regularizations are equivalent. 
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Fig 1. Phase diagram for LT— >0. jxi^ ore horizontal and vertical axes. Green lines 
Ji2 = ±A*i are boundaries of the fundamental domain. Solid blue and red curves are the 
lines of first order phase transitions. The phase between the red lines has a 2 7^ 0. The 
phase between the blue and the green lines has p 7^ 0. All other condensates vanish. 



thermodynamically preferred. An interesting limit of (2.24) is c— >oo, which corresponds 
to the large values of 

H(c) ^^ci (2.26) 

which corresponds to 

5S « H (2-27) 

7c 2 2 2 

Another possible configuration is the one with partially or completely restored chiral sym- 
metry. Consider the phase where o~\ = 0, 02 7^ 0. For this configuration we need to use 
(2.10) with djjX A = which corresponds to the limit c—>-oo. Hence, we can use 

some of the results above in the limit c— >oc. The value of chemical potential is now 

Hi = cici (2.28) 

which should be compared with (2.13) and (2.16). The value of SS is 

SSi r = El (g 7 (0)H(0) + (2-29) 

Analogously, the value of the action for the phase with completely restored chiral symmetry 
is 

SS 12r = E^(Jlj + J1I) (2.30) 



-120 



-100 



-20 



-40 



-60 



-80 




Fig 2. Phase diagram for LT— >0. p\^, are horizontal and vertical axes. Green line 
p,2 = — pi is a boundary of the fundamental domain. Blue curve is the line of first order 
phase transitions between the phase with p ^ (below) and p = (above). The leading 
asymptotics is p2 ~ — Pi- 

Comparing (2.29), (2.30) and (2.27) we deduce the asymptotic behavior of the line of the 
phase transitions between the phase with p 7^ and the phase with p = 0. This line (shown 
as a blue curve in Figs. 1,2) must behave asymptotically as JI2 ~ — pi at large values of 
Pi- 

More refined features of the phase diagram need numerical work. The results are sum- 
marized in Figs. 1,2. The green lines, P2 = ±/ii define the boundaries of the fundamental 
domain. For pi^ > p c the phase with straight branes dominates and all condensates van- 
ish. Note that p c is the value of ps at the point of the phase transition in the pi = 
plane. It is computed numerically to be p c = 0.11 at LT— >0. For — p c < p2 < p c but fx\ 
to the right of the blue line, the phase with straight brane "1" but curved brane "2" (with 
vanishing electric field) dominates. This is the phase with o~± =0, 02 7^ 0. Below the blue 
line is the habitat of the curved solution with p 7^ 0, o~i = 02 = 0. Finally, all condensates 
vanish between the blue and the red lines, for sufficiently large values of ps and pj. 

Fig. 2. depicts the phase diagram at a larger scale. As explained above, the blue 
phase transition line has the leading asymptotics P2 ~ —pi- The correction seems to 
behave like yfpl. 
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3. Intermediate temperatures 

In this Section we consider the case of finite LT. As we will see, for sufficiently large 
values of LT the phases with nonzero condensates do not exist. We also assume that the 
temperature is sufficiently high, so that gluons are deconfined. This means that the black 
hole has formed in the holographic dual. The metric is now 

ds 2 = 2 (f(U)dt 2 + SydXW + (dX^f) + (|) 2 + U 2 dQ\] (3.1) 

where 

m = l-§, U T= (J^1 (3 ,) 

The dilaton is given by (2.5). The analysis of the previous Section corresponds to the 
Ut = limit. To determine the phase diagram in the — ii<i plane at finite temperature 
we need to repeat the analysis of Section 2. 
The DBI action now reads 



S = 2 J dUUi ^ 1 - (duA )2 + f{U) 0Q {duX*)* (3.3) 

The analogs of (2.11) and (2.12) are 

(B X 4 )2 = U3(US + c 2 )R*f(Uo) 

[U ' U*{U* + <*)P{U)-Uji{US + <*)U*f{U )f{U) { ' } 

and 

m A) 2_ c 2 U 3 f(U) . 

[ ° uAo) U*(U* + U$*)f{U) - U*(l+V)f(U ) 1 J 

where we again used c 2 = UqC 2 . It will be convenient to introduce dimensionless variables 
by rescaling 

U Q = u E ; ^ = JiiEo (3.6) 



where E = 4tt\/L 2 as before. 
Another useful quantity is 



t = *f (3.7) 



Using equations of motion to compute L together with (3.6), and going to the rescaled 
variables, we obtain the following relation: 

t = y*G(c,y) (3.8) 
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Fig 3. t = y?G(c,y) [see eq. (3.8)] as a function of y at c = [red line] and c = 
1 x 10 6 [blue line]. There is a limiting temperature beyond which the phase with connected 
branes does not exists. 



where 



UtV t 2 



Uq J u 



where we used (3.2), (3.6) and (3.7), and 



X 



(i + c 2 )(i - y 2 ) 

G( C ,y)= I dx X \ (x 5 +S 2 )(a .3 _ ^3)2 _ (1 + c2)( x 3 _ ?/ 3)( 1 _ y 3) ( 3 - 10 ) 
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Fig 4. Phase transition line at fij = 0. Horizontal axis is \i c . Vertical axis is t. 
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Fig 5. Phase diagram fort = 0.2. ^1,2 are horizontal and vertical axes. Green lines 
JI2 = ±/5i are boundaries of the fundamental domain. Solid blue and red curves are the 
lines of first order phase transitions. The phase between the red lines has a 2 7^ 0. The 
phase between the blue and the green lines has p 7^ 0. All other condensates vanish. 



P2 




Fig 6. Phase diagram fort = 0.2. ore horizontal and vertical axes. Green line 
JI2 = —pi is a boundary of the fundamental domain. Blue curve is the line of first order 
phase transitions between the phase with p 7^ (below) and p = (above). The leading 
asymptotics is P2 ~ — A*i- 

At a given value of t (or, equivalently, LT) (3.8) determines y (and uq) as a function of c. 
The behavior of the right hand side of (3.8) is shown in Fig. 3. There are two solutions of 
eq. (3.8). The one with larger y comes closer to the black hole horizon than its smaller y 
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Fig 7. line of phase transition [shown in blue; bottom line] in the [xb = plane. 
Horizontal axis is 2/Zj; vertical axis is t. Red [top] line is the limiting value oft. 



counterpart. It is the latter one which connects to the vacuum solution at T = and it is 
not hard to check that it is thermodynamically preferred. In the following we will focus on 
this solution. Note that for any value of fii there is a limiting temperature beyond which 
the curved brane solution does not exists. 
The next step is to compute fxj 

1 ~ t 2 

d^i — M2I = — F(c,y) (3.11) 

2 y 



where 




F(c,y)= I dx x \ (a;5 +B2)(x s _ y3) _ (1 + c2)(l - ya) ( 3 - 12 ) 



Using (3.8) and (3.11) one can determine c and y in terms of Ji\ t 2- To determine the 
thermodynamically preferred state we need to compute the action in various phases. This 
is again done in a similar way as in the previous Section. The expression for SS with 
non-vanishing p is [compare with (2.22)] 



5S 



u T Juo 1 V U ^ U " + c2 )^) - ^o 5 + c 2 )f(U ) 



(3.13) 

In terms of rescaled variables, this becomes 

6S = 2E$—H(c,y) (3.14) 

13 



where 



2 7 f°° <- ( l~ x 5 (x 3 -v 3 ) \ 

H(c,y) = -(1 -y*)+ / dtf l-< - — — — - — i— — ;/ — (3.15) 

7 Ji \ V (x 5 +&)(x 3 -y 3 ) - (l+2 2 )(l-y 3 )J 

Note that the functions G(c, y), F(c, y), if(c, y) defined in (3.10), (3.12) and (3.15) respec- 
tively, reduce to the functions G(c), F(c), H(c), which were used in the previous Section, 
in the limit y^O. 

The value of SS in the phase with curved branes and zero electric field (<ti = 02 ^ 
0, p = 0) is 

r V 

5S = 2E*—H(0,y) (3.16) 

where y is determined from (3.8) with c = 0. To analyze the phase diagram we also need 
to compute the value of the action in the phases with Oi = 0, cr^ 7^ and <j\ ^ = 0. The 
former is given by 

SStr= ^ + jy U i(i- ] f^pj (3,7) 

where c is related to fa as 



z 2 

In terms of the rescaled quantities, 



^ = J^dU ] f^TT, (3.18) 



_*s , 2 r ( [~J r(A)r(f) c A.1.6. ?°\\ 

^-^- + ^ t ^1+^-1 + ca - ? = ^F^-,-,-,-^^ 

(3.19) 



where c is determined through 



^2 



00 



Finally, the value of the SS for the phase with vanishing condensates is 

SS\2r = (JS'lr- + 55 < 2r — 255'o (3-21) 

We are ready to analyze the phase diagram. As explained at the end of Section 2, J1 C (T) 
defines the line of phase transitions in the /Ib, T plane at Jlj = 0. Here we reproduce this 
plot for completeness in Fig. 4. 



14 



Figs. 5,6 contain the slice of the phase diagram in the Ji\ — J12 plane for a sample value 
of t = 0.2. The pictures look similar to the ones in Figs. 1,2. The only visible difference 
is shrinking of the Jx c with temperature, described by the curve in Fig. 4. Therefore the 
picture which emerges at large Jii is the surface of first order phase transitions in the 
JiB,Jii,T space whose sections are the blue lines in Figs. 2, 6. This surface resembles 
hyperboloid with Jii being the axis of rotation. To determine the behavior of the surface 
at large Jlj, we consider Jib = plane. 

In Fig. 7 we plot the curve which is obtained by intersecting the phase transition 
surface with the Jib = plane. Red line on this graph denotes the limiting value of t, 
beyond which the curved brane solution does not exist. Both lines seem to go to limiting 
values of t at large From this picture we see that the phase transition surface never 
intersects the limiting surface, and the transition remains first order. 

4. Conclusions 

In this paper we analyzed the phase diagram of holographic QCD at finite temperature 
and baryon and isospin chemical potentials, defined by (2.2). The action (2.1) is invariant 
under [i\ <-> //2, which corresponds to [ii — > — /zj. Moreover, the physics is also invariant 
under y^— > — /Zj. Hence, the fundamental domain in the \i<i plane can be taken to be 
(1.1). We found an intricate structure with various phases separated by the first order 
phase transitions. 5 . 

In particular, at large values of we found a surface of first order phase transitions 
which separates the phase with the nonzero value of p existing at smaller values of \ib 
from the chirally restored phase at larger values of \ib- 

It is not hard to reinstate gluon confinement /deconfinement transition in the phase 
diagram. It is represented as a surface T = l/2ivR±. For temperatures lower than this 
value chiral symmetry is necessarily broken and the system is in the phase with nonzero 
electric flux on the branes. We do not analyze this phase in detail here. It is clear that in 
the regions of the phase space where chiral symmetry is broken, the surface T = l/2itR± 
describes the coincident confinement/deconfinement and chiral phase transitions. 

5 As usual in the DBI analysis, the free energies of different phases have different values of the 
derivatives at the crossing point, indicating first order phase transitions. 
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An important observation concerns the energy scales in the problem. We have seen 
that the relevant variables in the phase space are LT and n/E . This is consistent with the 
energy scales of the mesons being 0(1/ L) and constituent masses of quarks being O(E ). 

The DBI action used to analyze the structure of the phase diagram in this paper is 
essentially abelian. It would be interesting to investigate the effect of non-abelian degrees of 
freedom, in particular with regards to the stability of the solutions. The existence of other 
phases also cannot be ruled out. Other directions for future research include generalizing 
our results to a higher number of flavors and studying other holographic models of QCD- 
like theories. 
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